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. The generic predictions of two-texture zero neutrino mass matrices in the flavor 

\^ I basis have been examined especiahy in relation to the degeneracies between mass ma- 

CD ' trices within a class and interesting constraints on the neutrino parameters have been 

obtained. It is shown that the knowledge of the octant of O23, the sign of cos 5 and 
neutrino mass hierarchy can be used to lift these degeneracies. 
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1 Introduction 



H ' Mass matrices provide important tools for the investigation of the underlying symmetries 
and the resulting dynamics. The first step in this direction is the reconstruction of the neu- 
trino mass matrix in the flavor basis. However, the reconstruction results in a large variety 
of possible structures of mass matrices depending strongly on the mass scale, mass hierarchy 
and the Majorana phases. However, the relatively weak dependence on some oscillation pa- 
rameters (^23 and 6) results in the degeneracy of possible neutrino mass matrices. The mass 
matrix for Majorana neutrinos contains nine physical parameters including the three mass 
eigenvalues, three mixing angles and the three CP-violating phases. The two squared-mass 
differences (Am^2 and Aml^) and the two mixing angles (6*12 and ^23) have been measured in 
solar, atmospheric and reactor experiments. The third mixing angle ^13 and the Dirac-type 
CP-violating phase S are expected to be measured in the forthcoming neutrino oscillation 
experiments. The possible measurement of the effective Majorana mass in neutrinoless dou- 
ble f3 decay searches will provide an additional constraint on the remaining three neutrino 
parameters viz. the neutrino mass scale and two Majorana-type CP-violating phases. While 
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the neutrino mass scale will be independently determined by the direct beta decay searches 
and cosmological observations, the two Majorana phases will not be uniquely determined 
from the measurement of effective Majorana mass even if the absolute neutrino mass scale 
is known. Under the circumstances, it is natural to employ other theoretical inputs for the 
reconstruction of the neutrino mass matrix. The possible form of these additional theoretical 
inputs are limited by the existing neutrino data. Several proposals have been made in the 
literature to restrict the form of the neutrino mass matrix and to reduce the number of free 
parameters which include presence of texture zeros [H [21 El IH [5] , requirement of zero deter- 
minant |6], the zero trace condition [7j to name just a few. However, the current neutrino 
oscillation data are consistent only with a limited number of texture schemes [U El [3l HI |5] . 
In particular, the current neutrino oscillation data disallow all neutrino mass matrices with 
three or more texture zeros in the flavor basis. Out of the fifteen possible neutrino mass 
matrices with two texture zeros, only seven are compatible with the current neutrino os- 
cillation data. The seven allowed two texture zero mass matrices have been classified into 
three categories. The two class A matrices of the types Ai and A2 give normal hierarchy 
(NH) of neutrino masses. The class B matrices of types Bi, B2, S3 and B4 yield a quasi- 
degenerate (QD) spectrum of neutrino masses. The single class C matrix corresponds to 
inverted hierarchy (IH) of neutrino masses. 

In the absence of a significant breakthrough in the theoretical understanding of the fermion 
flavors, the phenomenological approaches are bound to play a crucial role in interpreting new 
experimental data on quark and lepton mixing. These approaches are expected to provide 
useful hints towards unraveling the dynamics of fermion mass generation, CP violation and 
identification of possible underlying symmetries of the lepton flavors from which realistic 
models of lepton mass generation and flavor mixing could be, hopefully, constructed. 

Even though the grand uniflcation on its own does not shed any light on the flavor problem, 
the Grand Unifled Theories (GUTs) provide the optimal framework in which possible solu- 
tions to the flavor problem could be embedded. This is because the GUTs predict deflnite 
group theoretical relations between the fermion mass matrices. For this purpose, it is useful 
to flnd out possible leading order forms of the neutrino mass matrix in a basis in which the 
charged lepton mass matrix is diagonal. Such forms of neutrino mass matrix provide useful 
hints for model building which will eventually shed important light on the dynamics of lep- 
ton mass generation and flavor mixing. For example, a phenomenologically favored texture 
of quark mass matrix has been presented earlier [8]. In the spirit of quark-lepton similar- 
ity, the same texture has been prescribed for the charged lepton and Dirac neutrino mass 
matrices. The same texture for the right handed neutrino mass matrix in the see-saw mech- 
anism might follow from universal flavor symmetry hidden in a more fundamental theory of 
mass generation. Thus, the texture zeros in different positions of the neutrino mass matrix, 
in particular and fermion mass matrices, in general could be consequence of an underlying 
symmetry [9l [10]. Such universal textures of fermion mass matrices can, theoretically, be 
obtained in the context of GUTs based on SO(IO) [11]. Moreover, neutrino mass matrices 
with texture zeros have important implications for leptogenesis |12j . 

In the present work, we examine phenomenological implications of all possible neutrino 
mass matrices with two texture zeros. Neutrino mass matrices within a class were thought 
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to have identical phenomenological consequences [H [21 [3] leading to degeneracy of mass 
matrices within a class. Thus, there is a two-fold degeneracy in class A while the neutrino 
mass matrices of class B exhibit an eight-fold degeneracy since normal/inverted hierarchies 
are practically indistinguishable because of quasi-degenerate spectrum of neutrino masses in 
this class. We study this degeneracy in detail and discuss possible ways to lift this degeneracy. 
It is found that the deviation of atmospheric mixing from maximality and the quadrant of 
the Dirac-type CP-violating phase 6 can be used to distinguish the mass matrices within a 
class. We, also, note that the determination of hierarchy will have important implications 
for class B neutrino mass matrices. The prospects for the measurement of ^13 for class A 
neutrino mass matrices are quite optimistic since a definite lower bound on ^13 is obtained 
for this class. For neutrino mass matrices of class B, CP- violation will be near maximal 
if ^13 is large. Definite lower bounds on the effective Majorana mass M^e are obtained for 
neutrino mass matrices of class B and C. Class D of neutrino mass matrices is disallowed in 
our analysis. 



2 Neutrino mass matrix 

The neutrino mass matrix, M, can be parameterized in terms of three neutrino mass eigen- 
values (mi, m2, m-3), three neutrino mixing angles {O12, O23, O13) and one Dirac-type CP 
violating phase, 6. If neutrinos are Majorana particles then there are two additional CP 
violating phases a, f3 in the neutrino mixing matrix. The complex symmetric mass matrix 
M can be diagonalized by a complex unitary matrix V: 

M = VAlf'^W^ (1) 

where M^*°^ = Diag{mi,m2,ms}- The neutrino mixing matrix V [13] 

/ C12C13 S12C13 5136"*'^ \ 

V = UP=\ -S12C23 - Ci2S23Sl3e*'^ C12C23 - 5125235136*"^ S23C13 

V 512523 - Ci2C235i3e*'' -Ci2523 - 5i2C235i3e*'' C23C13 / 

where Sij = sin 6ij and Cij = cos 6ij. The matrix V is called the neutrino mixing matrix or 
PMNS matrix. The matrix U is the lepton analogue of the CKM quark mixing matrix and 
P contains the two Majorana phases. 

The elements of the neutrino mass matrix can be calculated from Eq. (1). Some of the 
elements of M, which can be equated to zero in the various allowed texture zero schemes, 
are given by 

Mee = cl^cl^mi + c?35?2"^2e^*" + sl^m^e^''^ , (3) 
= Ci3{5i3S23e^'^(e2^^m3 - sl^e^'''m2) - 012623512(^1 - e^^^ma) - c22Si3S23e^^mi}, (4) 
M^r = ci3{si3C23e'\e^''^m3 - 5?2e^'""^2) + Ci25235i2(mi - e^^^ma) - c?2Si3C23e*'^mi}, (5) 

= mi(c235i2 + e*^Ci25i3523)^ + e2*"m2 (C12C23 - 6^-^512513523)^ + e2*(^+'^) 77130^35^3 (6) 



can be written as 
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Type 




A, 


M — n A/f — n 


An 




Bi 


Me. = 0, M^^ = 


B2 


Me^ = 0, Mrr = 


B3 


= 0, M^^ = 


B, 


Me. = 0, M.. = 


C 


M^^ = 0, M.. = 


Di 


M^^ = 0, M^. = 


D2 


Mrr = 0, M^r = 



Table 1: Allowed two texture zero mass matrices. 

and 

Mrr = mi(s23Sl2 - Cu S 13023)^ + e^'^ms (Ci2S23 + e*^Si2Si3C23)^ + c^^^^^'^) rus^^ . (7) 

It will be helpful to note from Eqs. (4-7) that the transformation 

T : 623^^-623,6-^6 + 71 (8) 

transforms M^^ to -M^r and M^^ to M... Therefore, if Mg^ vanishes for ^23 and 6, then M^r 
vanishes for ^ — ^23 and 6 + 11. Similarly, if M^^ vanishes for ^23 and 6, then M.. vanishes 
for I — ^23 and 6 + tt. This fact can be, gainfully, used for distinguishing the subcategories 
of neutrino mass matrices within a class. 



3 Two texture zeros 

The neutrino mass matrix, M, in the charged lepton basis is given by Eq. (1). Out of the 
fifteen possible patterns for two texture zeros in M, the seven possibilities [HE], which are 
found to be consistent with the current neutrino oscillation data, are listed in Table 1. In 
addition, we have, also, included the two mass matrices of class D in the list. These mass 
matrices were disallowed in some of the earlier analyses [H [2]. However, it was found in a 
subsequent detailed numerical analysis by Guo and Xing [3] that the neutrino mass matrices 
of class D are, marginally, allowed. The expressions for the vanishing elements in different 
two- zero texture schemes have been listed in Eqs. (4-7). 

The transformation T defined in Eq. (8) transforms neutrino mass matrices of type Ai to 
A2 and neutrino mass matrices of type Bi {B4) to B2 (-B3). Therefore, the predictions for 
neutrino mass matrices of types Ai and A2 will be identical for all neutrino parameters except 
^23 and/or 6. Similarly, the predictions of neutrino mass matrices for all the four types of 
mass matrices in class B will be identical for all neutrino parameters except 6*23 and/or 6. 
However, the different types of mass matrices will have different predictions for the octant 
of 6*23 and the sign of cos 6. 
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The two texture zeros in the neutrino mass matrix give two complex equations viz. 



Mab = 0, = 



(9) 



where a, b, p and q can take values e, /i and r. Eqs. (9) can, also, be written as 

miUaiUn + m2^7„2t/(,2e'^" + m3Ua3Uh3e^'^^+'^ = (10) 



and 



(11) 



where U has been defined in Eq. (2). These two complex equations involve nine physical 
parameters mi, m2 , ma, ^12, ^23, ^13 and CP-violating phases a, and 5. The two mixing 
angles (6*12, 6*23) and two mass-squared difi^erences (Ami2, Am^a) are known from the solar, 
atmospheric and reactor neutrino experiments. The masses m2 and ma can be calculated 
from the mass-squared differences Ami2 and Am^g using the relations 



m2 = ^Jm^~-\-~Am 
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and 



ma = Y^mi + Am^a- 



(12) 



(13) 



Thus, we have two complex relations relating five unknown parameters viz. mi, ^13, a, f3 
and 5. Therefore, if one out of these five parameters is assumed, other four parameters can 
be predicted. Thus, neutrino mass matrices with two texture zeros in the fiavor basis have 

strong predictive power. 

Solving Eqs. (10) and (11) simultaneously, we obtain 



and 



ma 



UalUblUp2Uq2 — Ua2Ub2UplUqi 

Up2Uq2UazUbZ — Ua2Ub2UpsUq3 



1712 



(14) 



(15) 



UalUblUpzUqZ — UazUbzUpiUqi 

Using Eqs. (14) and (15), the two mass ratios ^) and Majorana phases (a, (5) can be 
written as 

Up2Uq2Ua3Ub3 — Ua2Ub2Up'iUqZ ^^^-^ 

(17) 
(18) 



mi 
m2 



mi\ 
ma/ 



UalUblUpsUqS 
UpzUq3Ua2Ub2 



UazUbzUplUqi 
UazUbzUp2Uq2 



a 



and 



—^arg 



1 



UalUblUp2Uq2 — t^tt2^62^9lt^yl 

/ Up2Uq2Ua3Ub3 " ^a2 f42 f^p3 ^53 
\UalUblUp3Uq3 — Ua3Ub3UpiUqi 

Up3Uq3Ua2Ub2 — Ua3Ub3Up2Uq2\ 



(3 — — arg 

2 \UalUblUp2Uq2 — Ua2Ub2UpiUqi 



-5. 



(19) 



5 



Eqs. (18) and (19) give the Majorana phases a and /3 in terms of ^13 and b since Q\2 and ^23 
are known experimentally. Similarly, Eqs. (16) and (17) give the mass ratios (^^, as 
functions of ^13 and b. Since, A'ml2 and Am^^ are known experimentally, the values of mass 
ratios ^) from Eqs. (16) and (17) can be used to calculate mi. This can be done by 
inverting Eqs. (12) and (13) to obtain the two values of mi, viz. 



mi 



mi\ 

1712/ 



\ 



Amfa 
1 _ 



(20) 



and 



mi 



mi\ 



\ 



Am?2 + Ami 



23 



1 - 



(21) 



The two values of mi obtained above from the mass ratios and respectively, 

must be equal. Thus, we can constrain (^13, 6) plane using the experimental inputs for 
Amfg, Am^g, 612 and 6123. The Dirac-type CP violating phase, 6, is given full variation from 
0° to 360° and 6*13 is varied from zero to its upper bound from CHOOZ experiment. We use 
the current best fit values of the oscillation parameters with 1, 2 and 3 a [H] errors given 
below: 



Ami2 
Am23 

■'12 



'23 



'13 



< 



7.9 X 10-^ eV' , 

±2.6 nl' X 10-=^ eV' 
0.30 
0.50 

0.012, 0.025, 0.040 



+0.2, 
-0.2, 
+0.02, 
-0.02, 

+0.06, 
-0.05, 



0.4, 0.6 
0.4, 0.6 ^ 
0.06, 0.10 
0.04, 0.06 5 
0.13, 0.18 
0.12, 0.16 ) 



(22) 



Here, ('— ') sign with the value of Am23 is for normal (inverted) hierarchy. The analysis 
|14j incorporates not only the latest long baseline data for Am^g from the MINOS collab- 
oration [15j but also the updated KamLAND and SNO data [16]. We vary the oscillation 
parameters within their 1, 2 and 3 a C.L. ranges with uniform distributions and calculate 
their predictions at various confidence levels. However, we report the correlations plots at 1 
a C.L.. We, also, calculate the Jarlskog rephasing invariant quantity [H] 



J = S12S23S13C12C23C13 sin 5 



(23) 



for the allowed parameter space. 

In our numerical analysis outlined above, the two squared-mass differences (Amf2 and Am^g) 
enter as two independent experimental inputs while the earlier analyses [H [3] require the 
ratio of two known mass-squared differences 



^ _ Am^2 
Ami3 



X _ (mi \2 



^ '"2 



(im )2 



(24) 
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Type 


mi 

my. 


mi 
ma 


a 




A, 


tan^ 612 


tan ^12 tan 6*23 sin 6*13 


n + ^TT 




A2 


tan^ 6^12 


tan 6^12 sin 6^13 
tan 023 




= (n + i)7r 


Bi 


1 


tan^ 023 


nir 


(3 + S = (n + '^)tt 


B2 


1 


1 

tan"^ 023 


rni 


(3 + S = {n + ^)7r 


B3 


1 


tan^ 6'23 


nvr 


P + S= {n + ^)7r 


B4 


1 


1 

tan"^ 02fi 


nvr 


(3 + S = {n + ^)7T 


C 


1 

tan-^ 01 


tan 01 7, tan 202.s sin 0i ;i 


n + ivr 




Di 


tan^ 01 9, 


tan y2;i 
tan 012 sin 013 


nn 


= (n + i)7r 


D2 


tan^ 01 


tan02fitan0i2sin013 


n+\'K 


/3 + 1 = 



Table 2: Zeroth order approximations for neutrino mass matrices with two texture zeros. 



to lie in the experimentally allowed range in order to constrain the other neutrino param- 
eters. Our analysis makes direct use of the two mass-squared differences and has more 
constraining power whereas the earlier analyses which make use of the ratio R^, lose some 
of the constraining power because this procedure does not use the full experimental infor- 
mation currently available with us in the shape of two mass-squared differences. Moreover, 
since is a function of mass ratios, it does not depend upon the absolute neutrino mass 
scale. It is obvious that the two mass ratios (^) and (^j may yield the experimentally 
allowed values of Ri, [Eq. (24)] for mutually inconsistent values of mi [Eqs. (20) and (21)]. 
Such mass ratios were allowed in the earlier analyses based upon Ri, [U El |3] . However, our 
analysis selects only those mass ratios for which the values of mi obtained from Eqs. (20) 
and (21) are identical and the ratio R^ for these values of the mass ratios is automatically 
restricted to lie in the allowed experimental range. Moreover, the definition of R^ 



Amfg |m2 — mil 

^23 



Ru — A 2 ~ I 2 21 (^^) 

Ami, \mi — mi\ 



used in many earlier analyses does not make use of the knowledge of solar mass hierarchy to 
constrain the neutrino parameter space. Our analysis disallows the class D of two-texture 
zero neutrino mass matrices which has been marginally allowed by Guo and Xing [3j. 



4 Results and discussion 

In this section, we shall present the results of our numerical analysis for neutrino mass 
matrices of classes A, B and C based upon the approach described in the previous section. 
However, it will be helpful to know the Taylor series expansion of Eqs. (16-19) in the powers 
of Si3 to appreciate the numerical results. The zeroth order terms of the Taylor series 
expansion have been tabulated in Table 2. 
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C.L. 




a(deg) 




?) 


6*13 (deg) 


J 




1 


2.6 - 4.5 


83 - 97 


-80 - 


80 


5.2 - 6.3 


-0.025 - 


0.025 


2 


2.1 - 8.9 


77 - 103 


-90 - 


90 


4.2 - 9.1 


-0.037 - 


0.037 


3 


1.7 - 13.4 


73 - 107 


-90 - 


90 


3.3 - 11.5 


-0.046 - 


0.046 



Table 3: The predictions for neutrino mass matrices of class A. 



1 a predictions 


5 


A, 

A2 


-150° - 150° 41.8° - 48.2° 
30° - 330° 41.8° - 48.2° 



Table 4: The predictions for 6 and 623 for neutrino mass matrices of types Ai and A2. 

4.1 Class A 

It can be seen from Table 2 that both the mass ratios — , — <^ 1 for neutrino mass matrices 
of class A. So, neutrino mass matrices of class A give hierarchical spectrum of neutrino 
masses. This case has, already, been analyzed analytically [5] and, here, we shall discuss the 
earlier results only for the sake of completeness. 

The results for mi, a, (3, 613 and J for class A mass matrices have been summarized in Table 
3 at various confidence levels. These quantities are the same for mass matrices of types Ai 
and A2. It can be seen from Table 3 that the 3 a lower bound on 6*13 is 3.3°. The range for the 
Majorana-type CP-violating phase /3 at 1 cr C.L. is found to be —80° - 80°. However, if the 
neutrino oscillation parameters are allowed to vary beyond their present 1.2 a C.L. ranges, 
the full range for P (—90° - 90°) is allowed. At one standard deviation, the allowed range of 
S is (—150° - 150°) for type Ai and (30° - 330°) for type A2 and the allowed range of 6*23 is 
(41.8° - 48.2°) for type Ai and A2 [Table 4]. Just like /3, no constraint on 6 is obtained above 

1.2 a C.L. The accuracies of the oscillation parameters required to distinguish between the 
theoretical predictions of neutrino mass matrices of types Ai and A2 depend on the upper 
bound on ^13. With the current precision of the oscillation parameters, the neutrino mass 
matrices of types Ai and A2 can be distinguished at 1 a, 2 a and 3 a C.L. for ^13 < 5.7°, 5.2° 
and 4.8° respectively. 

Fig. 1 (identical for Ai and A2) depicts the correlation plots of mi, a, P and J at one 
standard deviation. Fig. 1(a) shows a as a function of mi and Figl(b) shows as a function 
of a. The large spread in the (mi, a) plot is due to the errors in the neutrino oscillation 
parameters. However, there is a very strong correlation between a and P as indicated in 
Table 2. Such a correlation between the Majorana phases was noted earlier [3]. However, 
full ranges (—90° - 90°) are allowed for the two Majorana phases in that analysis [3]. In 
contrast, we obtain a very narrow range for the Majorana phase a around 90° [Cf. Table 3]. 
Similarly, the Majorana phase P is, also, constrained if the oscillation parameters are limited 
to their 1 a ranges. Fig. 1(c) depicts the correlation between J and mi. 

In Fig. 2, we depict the correlation plots of a and P with 6 for matrices of types Ai (left 
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panel) and A2 (right panel). We, also, show the correlation plots of S and 623 with one 
another as well as with ^13 in Fig. 2. The Dirac-type CP-violating phase 6 is, strongly, 
correlated with the Majorana-type CP- violating phase (3 [Cf. Table 2]. It can be seen from 
(a, (5) correlation plots [Fig. 2(a) and 2(b)] and (/3,5) correlation plots [Fig. 2(c) and 2(d)] 
that there are small deviations in the values of a around 90*^ and the correlation between P 
and 6 is, almost, linear. The fact that 6+2/3 ~ O'' for Ai type mass matrices and 6+2(3 ~ 180" 
for A2 type mass matrices [cf. Table 2] is apparent from the (/5,5) correlation plots [Fig. 2(c) 
and 2(d)]. The {(3,6) correlation was noted earlier by Xing [2] in its approximate form in a 
different parameterization. The {6, 623) plots, clearly, illustrate the point that the neutrino 
mass matrices of types Ai and A2 prefer different regions on the {6, 623) plane. This is 
contrary to the analysis done by Xing [3J where no constraints on 6 and 623 have been 
obtained. In fact, this feature is crucial for distinguishing mass matrices of types Ai and 
A2 which were found to be degenerate in the earlier analyses [H [21 [3] . The constraints on 6 
and 623 are very sensitive to the values of ^13. For the values of ^13 smaller than la CHOOZ 
bound, the constraints on 6 and 623 become stronger which can be seen from the (^13, ^23) and 
(6*13,5) correlation plots [Figures 2(g-j)]. For example, if ^^13 > 5.6°, then 6*23 > 46° (above 
maximal) for type Ai and ^23 < 44° (below maximal) for type A2. It can, also, be seen from 
the (^13, ^23) correlation plots in Fig. 2 that the deviation of ^23 from maximality is larger 
for smaller values of ^13. Therefore, if future experiments measure ^13 below its present la 
bound, the neutrino mass matrices of types Ai and A2 will have different predictions for 6 
and 6*23 with no overlap. However, it would be difficult to differentiate between matrices of 
types Ai and A2 if ^13 is found to be above its present la range. As noted earlier, different 
quadrants for 6 are selected for neutrino mass matrices of types Ai and A2. It can, also, be 
seen in Fig. 2 that the points 5 = 0, tt are ruled out implying that neutrino mass matrices of 
class A are necessarily CP-violating. This feature is, also, apparent in Fig. 2 which depict 
J as a function of 6 for types Ai [Fig. 2(k)] and A2 [Fig. 2(1)]. 



4.2 Class B 

It can be seen from Table 2 that the mass ratios — and — are approximately equal to 
one for neutrino mass matrices of class B. So, neutrino mass matrices of class B give quasi- 
degenerate spectrum of neutrino masses. For neutrino mass matrices of class B, the ratio 
upto first order in S13, is given by 

151 : — = 1H 0-S13COSO, 

"1-2 Ci2Si2S^3 
D . S23 . 

B2 : — = 1 3-Si3COSd, 

^^2 C12S12C23 

mi ^ C23tan^g23 

B3 : — = 1 3— Si3COSd, 

m2 C12S12S23 

mi g23 COt^ ^23 r 

Bi : — = 1H 3 — S13COSO. (z6j 

m2 C12S12C23 

Since, the mass ratio ^ is always smaller than one, we find that cos 5 should be negative 
for Bi and B^ and positive for B2 and S3, a result that can be gainfully used to distinguish 
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C.L. 


a 






J 


1 a 


-2.0'^- 


■2.0" 


> 0.064eV 


-0.025-0.025 


2 a 


-7.7°- 


-7.7° 


> 0.037eV" 


-0.036-0.036 


3 a 


-22°- 


-22° 


> 0.023eV" 


-0.045-0.045 



Table 5: The predictions for neutrino mass matrices of class B. 

between mass matrices of class B. Hence, Bi and (or B2 and ,83) will have identical 
predictions for 5. Similarly, it can be seen that Bi and B^ (or B2 and S4) will have, almost, 
identical predictions for 6*23 if it is near maximality. However, Bi and B^ (or B2 and B^) can 
be distinguished from each other by their predictions for the octant of ^23 and Bi and B^ 
(or B2 and B4) can be distinguished from each other by their predictions for the quadrant 
of 5. Thus, the four-fold degeneracy in class B is, now, lifted. Moreover, it follows from 
Eqs. (26) that the second term should be suppressed by S13 cos 5 if neutrino mass spectrum 
is quasi-degenerate (mi ~ 777,2). For large 6*13, the Dirac phase 5 should be peaked around 
90° or 270° to keep the term si3Cos5 small. However, for small ^13, no additional constraint 
on 5 is obtained. The fact that cos 5 should be negative for Bi and i?4 type mass matrices 
and positive for B2 and B^ type mass matrices for the correct sign of Aml2 results in the 
elimination of two quadrants of S for each of these subgroups irrespective of the values of the 
oscillation parameters and their errors. In fact, the earlier analyses [21 [3] do not constrain 
the Dirac phase S at all for these subgroup of mass matrices. In the analysis by Xing [2], 
the mass ratios — and — were examined and this important constraint on 6 coming from 

ma ma ^ ^ 

the mass ratio was comletely missed as a result. Another detailed numerical analysis 

by Guo and Xing [3j based upon R^, defined in Eq. (25) does not make use of the solar mass 
hierarchy as a result of which no constraints on 6 were obtained for neutrino mass matrices 
of class B. 

Now, we present the numerical results of our analysis for mass matrices of class B. The 
allowed range of 6 is (90° - 270°) for types Bi and B4 and (—90° - 90°) for types B2 and 
i?3. In other words, cos 5 is negative for types Bi and and positive for types B2 and 
i?3, as expected. This important result is valid irrespective of the ranges of the neutrino 
oscillation parameters entering in the analysis and is a consequence of the fact that the sign 
of Aml2 is positive. If 623 < 45°, Bi and B3 give normal hierarchy while B2 and B^ give 
inverted hierarchy. Similarly, if 623 > 45°, then Bi and B3 give inverted hierarchy while B2 
and i?4 give normal hierarchy. Again, this important result is true irrespective of the values 
of the neutrino oscillation parameters (and their errors) and follows from the dependence of 
mass ratio on tan^ ^23 [Cf. Table 2]. As noted earlier, the four types of neutrino mass 
matrices within class B differ in their predictions for 6 and ^23- There is no lower bound on 
613 in class B. The Majorana-type CP-violating phase /?, also, remains unconstrained. The 
predictions for a, Mee and J for class B mass matrices have been summarized in Table 5. 
These quantities are the same for mass matrices of types Bi, B2, B3 and B4. It can, also, 
be seen that a is restricted to a very small range around zero and definite lower bounds on 
Mee are obtained. It is important to note that 6*23 = 45° is disallowed for mass matrices of 
class B. 
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In Fig. 3, the correlation plots for all eight possible cases of B type mass matrices {Bi, 
B2, -B3 and -B4 with normal as well as inverted hierarchy) have been shown for a, P, 6, 613, 
mi, Mee and J. These plots show that B type mass matrices exhibit eightfold degeneracy 
in the sense that eight different types of mass matrices have identical predictions for these 
parameters. Figures 3(a-d) depict the correlations of a, (3, J and Mee with ^13. It can be 
seen from Fig. 3(a) that an upper bound on 6*13 constrains a to a very narrow range. Even at 
three standard deviations, the allowed range of a is —22° < a < 22° [Cf. Table 5]. However, 
no constraints are obtained on the Majorana phase f3 [Fig. 3(b)]. The correlation plot of J 
with 6'i3 shows a linear behaviour which is a consequence of ^13 being small and 6 being near 
90° or 270° for most of the allowed neutrino parameter space. Since, Mee — i^i, we give the 
correlation plots for Mee instead of mi in Fig. 3(d) where it can be seen that the absolute 
neutrino mass scale and 1-3 mixing angles are anti-correlated with each other. Next, we 
depict the correlations of the quantities Mee and J with the Majorana phases a and (3 in 
figures 3(e-h). The effective Majorana mass Mee diverges at a = and /? = and so these 
points are not allowed. A lower bound on Mee is required to constrain the Majorana phase 
/3 which, otherwise, remains unconstrained for mass matrices of class B [Fig. 3(f)]. It can 
be seen from Fig. 3(g) that J and a are directly correlated with each other. This is due to 
the fact that a is constrained to a very narrow range around zero degree so that sin a ~ a. 
The correlation between the two Majorana phases a and p is depicted in Fig. 3(i) which 
shows that (3 will be nearly zero for large a and vice- versa. Fig. 3(j) depicts the correlation 
between J and Mee- One can see that maximal CP- violation is possible near the lowest 
possible value of Mee in class B. It may be worthwhile to mention that the earlier analysis 
[3] fails to obtain any constraints on the ranges of CP-violating phases {a, (3 and 5) for the 
reasons discussed earher. 

Fig. 4 depicts the correlations which lift the degeneracy between the neutrino mass matrices 
of types Bi (S4) and B2 (B^). Since, the range of Dirac-type CP- violating phase, 6, is 
different for Bi {B4) and B2 (-B3), the eightfold degeneracy in B type mass matrices has, 
now, been reduced to four-fold. The left panel depicts the correlation plots for both Bi 
and B4 and the right panel depicts the correlation plots for B2 and B^. It can be seen 
from Fig. 4(a,b) that the range of Dirac-type CP-violating phase, 6, is 90° < 6 < 270° 
(—90° < <5 < 90°) for Bi and {B2 and -B3) type mass matrices, as expected. Also, 
there is a strong correlation between Majorana-type CP violating phase, f3, and Dirac-type 
CP-violating phase, 6 reinforcing the zeroth order result (^13 + 6 = (^+1)^1^ given in Table 
2. It can be seen from Fig. 4(c,d) that ^13 remains unconstrained in all cases because full 
range of ^13 is allowed at 5 = 90° and 270°. However, as we deviate from these values of 
6, 613 decreases rapidly to very small values. So, if the CP violation is found to be non- 
maximal, ^13 will be constrained to very small values. However, if CP violation is found to 
be nearly maximal, ^13 can be large. Thus, the suppression factor S13 cos 6 in the first order 
correction term in the Taylor expansion for is small as expected from the zeroth order 
approximation for this class. 

Fig. 5 depicts the correlations of Mee and a with 623. The correlation plots with 623 can 
be used to further reduce the remaining four-fold degeneracy in B type mass matrices to a 
two-fold degeneracy which will be lifted by the determination of the neutrino mass hierarchy. 
The left panel (right panel) depicts the correlation plots for Bi{NH) and B2{IH) {Bi{IH) 
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(3 J M, 



1 a 

2 (J 

3 a 



-14"-14" -0.024-0.024 > 0.016 eV 
-16°-16° -0.034-0.034 > 0.013 eV 
-18°-18° -0.045-0.045 > 0.009 eV 



Table 6: The predictions for /3, J and Mee for neutrino mass matrices of class C. 

and B2{NH)). The correlation plots of M^e and a with 6*23 for B^{NH) and B^^IH) (or 
Bs{IH) and Bi{NH)) will be, almost, identical to the plots given in Fig. 5 because Bi 
and i?3 (or B2 and ^4) have, almost, identical predictions for ^23 ^-nd, therefore, are not 
given here. It can be seen from Fig. 5(a,b) that the effective Majorana mass Mgg is strongly 
correlated with 62^. Moreover, maximal value of ^^23 is not allowed by the current neutrino 
oscillation data since M^e diverges at 6*23 = f. Larger the deviation of ^23 from maximality, 
smaller is the value of M^e- Also, a large value of a imphes large deviations of ^23 from 
maximality. The extent of the deviation of 2-3 mixing from maximality is governed by the 
magnitude of the upper bound on M^e while the direction of the deviation from maximality 
is governed by the quadrant of 5 and the neutrino mass hierarchy. 



It can be seen from Table 2 that the mass ratio — is greater than one at zeroth order. 



However, the higher order terms in the Taylor series expansion of the ratio ^ may make 
it smaller than one. Detailed numerical analysis shows that this actually happens and the 
mass matrices of class C are marginally allowed. The Taylor series expansion for ^ to first 
order in S13 is given by 



For ^ to be smaller than one, the term (tan 2^23 cos 5) should be positive. Therefore, 
when ^23 < 45°, we must have —90° < 5 < 90°. Similarly, when ^23 > 45°, we must 
have 90° < 5 < 270°. Just like class B, the correlation between the octant of ^23 and 
the quadrant of 5 obtained here for class C is independent of the ranges of the neutrino 
oscillation parameters and is a generic feature of two texture zero mass matrices. The points 
6 = 90°, 270° and ^13 = 0° are not allowed because the first order correction term vanishes 
at these points. The mixing angle ^23 should approach 45° as 5 approaches 90° or 270° and 
^13 approaches 0° so that the term (si3 tan 2^23 cos S) can make the ratio ^ smaller than 
one. The mass ratio — is greater than one at the zeroth order. Therefore, neutrino mass 

m3 o ' 

matrices of type C may be expected to be consistent with inverted hierarchy only. However, 
we shall see in the detailed numerical analysis that the normal hierarchy is, also, allowed 
marginally. 

For mass matrices of type C, with inverted hierarchy, the allowed ranges for the parameters P, 
J and ^-re given in Table 6. All other parameters remain unconstrained. The correlation 



4.3 Class C 




(27) 
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plots between various parameters have been depicted in Figs. 6-8. In Fig. 6, we have given 
the correlations of a, P, ^13, 6*23, M^e and J with 6. As discussed earlier, the values S = 90° 
and 270*^ are not allowed [Fig. 6(a)]. Therefore, maximal CP- violation is not allowed for 
mass matrices of class C. For the maximal values of 6, a should be equal to zero. Therefore, 
the point a = 0*^ is, also, not allowed. However, all other values of a are allowed. It can 
be seen from Fig. 6(b) that the Majorana phase /3 is a periodic function of 6 with a period 
of 180° and 13 can be zero at 5 = 0° or 180°. However, the points (/5 = 0°, 5 = 90°) and 
(/5 = 0°, 5 = 270°) are disallowed. The correlation between S and 613 has been depicted in 
Fig. 6(c). The correlation between S and 623 has been depicted in Fig. 6(d) where it can 
be seen that the maximal 2-3 mixing is not allowed. However, the mixing angle 623 will be 
below maximality if cos 5 is positive (— 90°-90°) and above maximality if cos 5 is negative 
(90°-270°). These conclusions are, completely, in agreement with the results obtained from 
the Taylor series expansion. Fig. 6(e) depicts the correlation between Mee and 5 which yields 
a lower bound of about 0.02 eV on M^e for 5 = 0° or 180° for neutrino mass matrices of 
class C. As S deviates from the values 0° and 180°, Mee increases and diverges as 6 becomes 
maximal. The correlation between J and 6 has been shown in Fig. 6(f). 

In Fig. 7, we have shown the correlations of a, /?, 623 and J with 6*13. In Fig. 7(a), (6*13, a) 
correlation has been depicted. Full range (— 90°-90°) for a is allowed for ^13 = 0. However, 
it can be seen from Fig. 7(a) that a lower bound on ^13 of about 4° will constrain a to the 
range (— 65°-65°). As noted earlier, f3 is, already, constrained to a small range around zero. 
However, it can be seen from Fig. 7(b) that P can be further constrained if ^13 is found to be 
larger than 3°. A lower bound on ^13 larger than 3° will disallow a finite range of (3 around 
zero. From the (6*13, ^23) correlation plot [Fig. 7(c)], it can be seen that larger the value 
of 6*13, larger will be the deviation of 623 from maximality. The point (^13 = 0°, 623 = 45°) 
is not allowed. The disallowed range of 623 around 45° increases with increase in ^13. For 
example, the disallowed range of 623 at 613 = 6° is approximately 44° — 46°. For small values 
of 6*13, J is proportional to ^13 and the constant of proportionality is determined by sin 5 [Cf. 
Eq. (23) ]. These features can be seen graphically in the correlation plot between J and 6*13 
depicted in Fig. 7(d). 

Some other important correlation plots have been shown in Fig. 8. It can be seen from the [a, 
(3) correlation plot [Fig. 8(a)] that the two Majorana phases can not vanish simultaneously. 
The correlation between Mge and a has been shown in Fig. 8(c) where it can be seen that 
an upper bound on M^e will reject a considerable range of a. The (Mge, (3) correlation plot 
has been given in Fig. 8(d). It can be seen that the allowed range of [3 (which is — 15°-15° 
for small values of Mee) can be further constrained if Mee is found to be larger than 0.03 eV. 
Fig. 8(e) depicts the correlation between Mee and ^23- It can be seen that Mee diverges for 
maximal 2-3 mixing. A lower bound on Mee can be used to constrain ^23- 

For C type mass matrices, with normal hierarchy Dirac-type CP-violation phase, 5, can 
have only two values 5 = 90° or 180°. Majorana-type CP-violating phases {a and f3) and 
^13 remain unconstrained. There is a strong correlation between ^13 and ^23 in this case. 
Larger the allowed range of ^13, larger will be the deviation of ^23 from maximality. The 
neutrino mass matrices of class C were examined analytically in [TD] where some of the 
results reported here for class C have been obtained. 
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4.4 Eightfold degeneracy 

The forthcoming neutrino experiments will aim at measuring the neutrino mass hierarchy, 
Dirac type CP violating phase 6 and the deviations of the atmospheric mixing angle from 
maximality [19]. The results of these experiments will fall in one of the following eight 
categories of phenomenological interest: 

1. NH, 6*23 < 45° and 90° < 5 < 270°, 

2. NH, 6*23 > 45° and -90° <6 < 90°, 

3. NH, 6'23 < 45° and -90° <6 < 90°, 

4. NH, ^23 > 45° and 90° < 5 < 270°, 

5. IH, ^23 > 45° and 90° < 5 < 270°, 

6. IH, ^23 < 45° and -90° <S < 90°, 

7. IH, ^23 > 45° and -90° <6 < 90°, 

8. IH, ^23 < 45° and 90° < 5 < 270°. 

The eight cases given above are degenerate for the present neutrino oscillation data because 
of the octant degeneracy of ^23 (i-e. if 6*23 < 45° or 6*23 > 45°), the intrinsic degeneracy in 
the sign of cos S and the two-fold degeneracy in the sign of Arn^^ . This degeneracy known as 
the eightfold degeneracy in the neutrino parameter space has been studied extensively in the 
literature fIU\ UTl [2^ . A specific project named Tokai-to-Kamioka-Korea (T2KK) two 
detector complex which will receive neutrino superbeams from J-PARC has been proposed 
[22] to resolve these degeneracies. The intrinsic degeneracy in the sign of cos 5 will be 
resolved by the spectrum information at T2KK and the degeneracy in the sign of Am^^ will 
be resolved by observing the difference in the earth matter effects between the intermediate 
and far detectors. The 6*23 octant degeneracy will be resolved by observing the difference in 
Am^2 oscillation effects between the intermediate and far detectors at T2KK. In addition, 
neutrino factories |23j and beta-beam experiments [24j have the potential to resolve the 
eightfold degeneracy. 

5 Conclusions 

If neutrinoless double beta decay searches give positive results and Mee is measured ex- 
perimentally or the atmospheric/reactor neutrino oscillation experiments confirm inverted 
hierarchy, the neutrino mass matrices of class A will be ruled out. A generic prediction of 
this class of models is a lower bound on 6*13 [Table 3]. If the forthcoming neutrino oscillation 
experiments measure ^13 below 3.3°, the neutrino mass matrices of class A will, again, be 
ruled out. If Mee is found to be non-zero, the mass matrices of class B or C may be allowed. 
The eight different types of mass matrices in class B can be distinguished from one another 
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Table 7: The consequences of resolving the eightfold experimental degeneracies for the neu- 
trino mass matrices with two texture zeros. 

by the future experiential data by the resolution of the eightfold degeneracy of the neutrino 
parameter space. The observation of correlations between various neutrino parameters like 
^13 5 ^23; ^ and Mee etc. will confirm the neutrino mass matrices with two texture zeros. 
Moreover, the CP-violation in lepton number conserving and non- conserving processes will 
be correlated in a definite way and there will be only one independent CP-violating measure 
for both types of processes. 

Main conclusions of the present study have been summarized in Table 7 which shows how 
the resolution of the degeneracies of neutrino parameter space will resolve the degeneracies 
of the neutrino mass matrices. As an illustration for the use of Table 7, we look at the results 
obtained for category 1 for which the signs of Amgs and cos ^23 are positive and the sign 
of cos 5 is negative. For this category, the neutrino mass matrices of type A2 and Bi (for 
normal hierarchy only) may explain the neutrino data whereas the mass matrices of class 
C are ruled out. Other possibilities summarized in Table 7 may be understood similarly. 
It is found that the degeneracies in the neutrino parameter space are inextricably linked 
to the degeneracies in the neutrino mass matrix and by finding the hierarchy, octant of 62^ 
and sign of cos 5, the future experiments will settle the fate of neutrino mass matrices with 
two texture zeros. The resolution of the eightfold degeneracy in the neutrino parameter 
space will, also, resolve the eightfold degeneracy of the neutrino mass matrices of class B 
as can be seen from figures 4 and 5. Figures 4 and 5, also, illustrate the decoupling of 
degeneracies in the sense that two mass matrices can be degenerate in their 5 predictions 
but may have different predictions for ^23 or vice- versa. A similar decoupling of degeneracies 
occur in neutrino oscillation experiments with matter effects [21] where an experiment may 
be insensitive to the sign of cos 5 but can resolve the degeneracy in the sign of cos 2^23 01 it 
can be insensitive to the sign of cos 2^23 but can resolve the degeneracy in the sign of cos 5. 
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Figure 1: Correlation plots for neutrino mass matrices of class A at one standard deviation. 
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Figure 2: Correlation plots for neutrino mass matrices of type Ai (left panel) and A2 (right 
panel). 
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Figure 4: The left panel shows the correlation plots for Bi or B4 while right panel shows 
the correlation plots for B2 or S3 type mass matrices. 
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Figure 6: The correlation plots for CP-violating phase S for neutrino mass matrices of type 
C with inverted hierarchy. 
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Figure 7: The correlation plots for the mixing angle ^13 for neutrino mass matrices of type 
C with inverted hierarchy. 
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Figure 8: Some other important correlation plots for neutrino mass matrices of type C with 
inverted hierarchy. 
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